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Abstract. This paper addresses the energy accumulation problem, in terms of the H2 norm, of 
linearly coupled dynamical networks. An interesting outer-coupling relationship is constructed, under 
■ which the H2 norm of the newly constructed network with column-input and row-output shaped 
matrices increases exponentially fast with the node number N: it increases generally much faster than 

On ■ 

2 when N is large while the H2 norm of each node is 1. However, the H2 norm of the network with 



a diffusive coupling is equal to 72-/V, i.e., increasing linearly, when the network is stable, where 72 is 
the H2 norm of a single node. And the H2 norm of the network with antisymmetrical coupling also 
increases, but rather slowly, with the node number N. Other networks with block-diagonal-input and 
block-diagonal-output matrices behave similarly. It demonstrates that the changes of H2 norms in 

5-H ' 



different networks are very complicated, despite the fact that the networks are linear. Finally, the 
influence of the H2 norm of the locally linearized network on the output of a network with Lur'e nodes 
is discussed. 
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1 Introduction 

Complex networks have attracted increasing attention from physicists, biologists, social scientists and 
control engineers in recent years [1-4, 9, 11-13]. Many complicated problems in various networks, 
^his work is supported by the National Science Foundation of China under grants 60674093, 60334030, 
10472001. 
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such as topological structures, small-world and scale-free characteristics, robustness and fragility, and 
self-similarity, have been studied (see [1, 4, 11-13] and references therein). In general, a complex 
dynamical network can be considered as a large-scale system with special interconnections among its 
dynamical nodes. The large-scale system theory has been extensively studied in the last three decades, 
and many interesting results have been established, on such basic issues as decentrally fixed modes, 
decentralized controllers design, diagonal Lyapunov function method, and M-matrix method |llj . 
Many complicated technical problems in large-scale systems, such as decomposition and aggregation, 
connective stability and instability, competitive equilibria, local instability, have also been studied (see 
E E] and references therein). New concepts such as 'small- world', 'scale- free', 'power law', etc. 
in the current studies of complex networks will in turn provide new opportunities to the traditional 
large-scale system theory. On the other hand, the mature large-scale system methods can provide 
effective tools to network stability analysis and control problems. Note that many dynamical aspects 
of networks such as network synchronization, network control, and dynamical modeling have been 
studied (see [2j O [10\ I13j). From a control theoretic viewpoint, a dynamical system has inputs and 
outputs, so does a network, where the input and output relation has to be considered. This paper 
attempts to explore the input to output H2 gain (i.e., H2 norm, which represents the energy of the 
system output with fixed input) of a general linear dynamical network with different couplings. The 
results show that the i^-norm changes can be unexpectively complicated, despite the linear structure 
of the network. 

The rest of this paper is organized as follows. In section 2, some preliminaries are introduced. In 
section 3, an interesting network with column-input and row-output matrices is constructed, whose 
energy increases exponentially fast with the node number N. In section 4, diffusively and antisym- 
metrically coupled networks are studied. Compared with the network studied in section 3, the H2 
norms of these network increases linearly and in particular the norm of the diffusive network is equal 
to the node number N. In section 5, networks with block-diagonal-input and block-diagonal-output 
matrices are similarly discussed. In section 6, for a Lur'e network, the influence of the energy of its 
locally linearized network on the original network output is studied. The last section concludes the 
paper. 

2 Preliminaries 

Consider a continuous-time linear system, 



£ = A x t,+B x w u yx = C& 
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where £ is the state of the system, w\ and y\ are the input and output of the system, respectively, 
and A\ G R nxn ,i?i G R nxm ,Ci G R' xn are given constant matrices. The transfer function from w\ 
to yi is G\(s) = C\{sl — A^^Bi. If A\ is stable, the H2 norm of system ([1]) is represented by the 
H2 norm of the transfer function Gi(s), which is defined by 



IG1OOII2 = \ ^~ I t^{GH{jw)Gx{jw)}dw. 



Y 

The H2 norm represents the power of the system output for fixed input, or the induced system gain, 
in control theory [16]. For computing the H2 norm, the following formula is convenient |16j . 
Lemma 1 If A\ is stable, then the H2 norm of system ([!]) is given by 



\\Gi(s)\\ 2 = ^tx a ce(BfQB 1 ), 
where Q satisfies the Lyapunov equation 

QAt + AjQ + Cfd = 0. 

□ 

When A\ is unstable, the L2 norm can be computed [16j . With a linear coupling, the N nodes, 
each described by ([1]), constitute a dynamical network as follows: 

N 

x i = A 1 x i + ^2 1ij A i2 x 3i i = l r--,N, (2) 

j = 1 

where 7^ G R, A\ is given as in ([T]), and A12 G R nxn is the inner coupling matrix describing the 
interconnections among components of Xj,j = 1,...,N. Let T = (^fij)NxN, where 7™ are given in 
([2]), which is referred to as the outer coupling matrix. Using the Kronecker product notation [9l IIP], 
a network with column-input and row-output shaped matrices can be rewritten as 

x={I N (g)A 1 + r(g)A 12 )x + Bw 1 , y = Cx, (3) 

where x = (xj, ■■■ , x T N ) T , B = E N (g) B u C = E% (g) C 1 and E N = (!,■■■ , 1) T G H N . 

To discuss the H2 norm of network ([3]), the following simple result for stability analysis is needed. 

Lemma 2 Suppose Aj, % = 1, ■ ■ ■ , k, are distinct eigenvalues of T, with k < N, not to count the 
multiplicity. Then, Ijy (g) A\ + T (g) A12 is stable if, and only if, A\ + \iA\2-, i = 1, • • • , k, are stable 
simultaneously. 

Proof Let T be a nonsingular matrix such that T _1 TT = J, where J is the Jordan form of T. 
Then, the following similarity transformation completes the proof easily: 

(r- 1 6?) i)(i N 6?) a 1 + r 6d a 12 )(t (&i) = i N 6§A 1 + j(9) a 12 . 



Throughout this paper, let Ai 2 
output inner coupling. 



□ 

= BiC\, which simply means that the network has an input and 



3 ii/2-norm energy accumulation 

In this section, consider the changes of H2 norm with the network size for a specially constructed 
network. For convenience, first consider the following N x N coupling matrix: 



r = {lij)NxN, la = 0; Jij = -1, if i > j, j < i - 2 and j is odd; otherwise, 7^ = 1. 
One can also define T inductively, that is, 



(4) 



N > 2, 



(5) 



r (Ar_i)x(iV-l) 71 
72 

where 71 is a column vector with all elements being 1, 72 is a row vector with the i-th element being 
-1 if i < N — 2 is odd and the other elements being 1. For example, when N = 5, one has 

/ 1 1 1 1 \ 
10 111 

r 5X 5 = -11011 
-11 1 01 
v -1 1 -1 1 j 

Lemma 3 When N is even, the distinct eigenvalues of T are -1 and 1, and T is similar to 



Aw 







-In 

2 

Jn 

2 



where In is the identity matrix of order y x ^ and Jn is a Jordan matrix of order y x y as 



J N_ 

2 



V 1 / 

When TV is odd, the distinct eigenvalues of T^vxTV are -1, 1 and 0, and Tjvxiv is similar to diag(Ajv-i, 0). 

Proof It suffices to prove this lemma by inductively determining the eigenvectors and generalized 
eigenvectors of T. 



4 



First, suppose that N is even. If N = 2, obviously, the eigenvalues of T 2x2 are -1 an d 1. If N = 4, 
one has 

-P 2 and T 4x4: Q 2 = Q2J2, 
T 

-1 1 \ / 1 1 




where P2 = \ I , Q2 — I I an d J2 is a Jordan matrix of order 

1 -1 J ' V — 1 —1 1 1 

2 as defined above. This shows that the statement holds for N = 4. Let Pn-2 and Pn (N > 6) be 

2 2 

N-2 „„j A7 ^ AT 



(iV — 2) x — ^- and iV x j matrices composing of column eigenvectors of T(n-2)x(n-2) an d ^NxN 
corresponding to the eigenvalue -1, respectively. Then, Pn_ is composed of (0, —1, 0, 1, lx (jv-4)) T 
and columns (O jv-2 x2 , Pn-2 ) T ■ 

Similarly, for N > 6, with Q2 as given above, let Q n-2 and Qn be (N — 2) x and N x 



matrices composing of column and generalized column eigenvectors of T(n-2)x(N-2) an d ^NxN 



N 
2 

Q N-2 

corresponding to the eigenvalue 1, respectively. Then, Qn is composed of I - I and gjv 



(—2a, —2a, f3\, ■ ■ ■ , /3n-4, 1> 1) T 5 where a is the first element of the last column of Q n-2 and in 



°2x^ 



2 



2 



* = !,••• , — 4, are determined by 



(rTVxTV — Inxn)Qn = 2(Q]v-2 , 0, 0) , 

2 ~^2~ 

where q jv-2 is the last column of Qn-2. 

2 2 

When ./V is odd, except a zero eigenvalue, it is similar to the case N being even. □ 

Remark 1 Column eigenvectors for the eigenvalue -1 and generalized column eigenvectors for the 
eigenvalue 1 are discussed in the proof of Lemma 3. Row and generalized row eigenvectors for the 
eigenvalues -1 and 1 can be studied similarly. 

With the coupled matrix as discussed above, the stability of network ([3]) can be easily tested. 

Corollary 1 If the number of nodes, N, is even, then network ([3]) is stable, i.e., A = I A\ + 
T (g) Ayi is stable if, and only if, A\ - A u and Ax + A X2 axe stable. If N is odd, A = I (g) Ai + T (g) A\ 2 
is stable if, and only if, A±, Ai — A\ 2 and A\ + A\ 2 are stable. □ 

Theorem 1 If the number of nodes, N, is even, the H 2 norm of system ([3]) is equal to the H 2 
norm of the following system: 

1 = {I®A 1 + J„®A 12 )Z + B„w i v = CnC, (6) 

where Jn is given as in Lemma 1, and Cn and Bn are determined inductively as follows: 

2 22 

Cn = (Cn,, (-nf-^f^Ci), with C 2 = {2d, 0), 

2 2 



in which, by denoting Bn__-, = (2BT , x\BT , x 2 BT , • • • , xn__„B7 , BT) T , 

2 2 



Bn = {2Bf, (2 + x 1 )Bj, (xi + x 2 )Bf, ■■■ , (x K „ + 1)B{ , Bf) T , with B 2 = (2Bf , Bjf. 

2 2 



If N(N > 5) is odd, the H 2 norm of system is equal to the H 2 norm of the following system: 

£ = Aoi + Bn±i q w, y = C t, (7) 

JV — 1 ~ 

where A Q = diag(/(g)Ai + J N -i (g) A 12 , 0), C a = (C N-i , (—1) 2 Ci), and B n+i is determined 

2 2 2 

inductively as follows: B n-i = (2B± , xiBf , x 2 Bf, • • • , x N-i_ 2 Bf , 2Bf , Bf) T , and 

S«±i = (2Bf , (2 + , (an + x 2 )Bf , • • • , (x«+i_ 2 + 2)Bf , 2Bf, Bf) T , with B 2o = (2Bf , 25f , Bff . 



Proof Taking a similarity transformation to change T into a Jordan form as in Lemma 3, one 
can complete the proof easily. □ 

Remark 2 From the Jordan matrix Jjv, one can see that the number 2 appears in it. This 

2 

means that the term 2 N will appear in the transfer function of system (|6j) or (|7|). It is imaginable 
that the H 2 norm of network ([3]) can be possibly very large. Of course, one may also take other 

similarity transformations such that the term 2 becomes 1 in Jn. However, it should be noted that 

2 

any similarity transformation does not change the system norm. 

Example 1 Consider the single node system (pQ) with A\ = —4.5, B\ = 1, C\ = 3. See Table 1 
for the H 2 norms of network J3]). 

Example 2 Consider the single node system ([TJ with 





f 1 \ 
















{ -4 -2.5 1 







, Ci 



2 2 



See Table 1 for the H 2 norms of network ([3]). 

Table 1 H2 norms of network ([3]) with different N. 



N 


1 


2 


3 


4 


5 




10 




20 




Ex 1 


1 


3.4641 


6.0828 


15.4919 


25.4755 




1.4363 x 10 3 




3.6864 x 10 6 




Ex 2 


1 


4.8990 


9.1761 


40.9878 


74.9765 




3.6676 x 10 4 




4.2018 x 10 9 





From Table 1, one can see that the H 2 norm of network ([3]) increases exponentially fast, much 
faster than 2 20 (= 1048576) when N is large. This means that even a single node system can generate 
a complicated network through simple coupling relations. 



4 Diffusive and antisymmetrical coupling relationships 

First, consider diffusive coupling networks. Change the outer coupling matrix r in as follows: 

Td = (lij)NxN, if i / j, 7ij are S iven as in © and7ii = 



A' 



^ 7y , i,j = 1,-.. ,N. (8) 



Obviously, the sum of every row of is zero, which generally refers to as diffusive coupling. 

Lemma 4 If N is odd, the eigenvalues of Tjj, given as in (JSj) , are —N,—(N — 2),—(N — 
2), • • • , —3, —3, —1, 0; especially, if N = 3, its eigenvalues are -3, -1 and 0. If N is even, the eigenvalues 
of T D are —N, —(N — 2), —(N — 2), • • • , -2, -2, 0; especially, if N = 2, its eigenvalues are -2 and 0. □ 

Corollary 2 Given A% and A12 with compatible dimensions, Ijy&lAi + Td(£)Ai2 is stable if, 
and only if, A±, A\ — A12, A\ — 3A\2, ■ ■ ■ , A\ — (N — 2)A±2 and A\ — NA\2 are all stable when N is 
odd, or Ai, Ai - 2A 12 , ■ ■ ■ , Ax — (N - 2)A 12 and A x - NA 12 are all stable when N is even. □ 

Theorem 2 If network ([3]) with r = Tjj is stable, then the H2 norm of ([3|) is 72 N, where 72 is 
the H2 norm of the single- node system (fTJ . 

Proof Because of the diffusive characteristic of Tp, is one of its eigenvalues and the corresponding 
eigenvector is In = (lj ■ • ■ , Let P be a nonsingular matrix such that P~ 1 TdP = J, where J is 
the Jordan form of Tp, and the last column of P is l^v and the last row of J is a zero vector. Let 
rj = (P^ 1 (g) I)x. Then, systems Q becomes 

V = (I N (g)A 1 + Jf^Au^ + iP-'lN^B^w, y = (l T N P(g)C 1 ) V . (9) 

Obviously, P -1 !^ = (0, • • • , 0, a) T , Let the last element of 1%P be (3. Then, one has 

a/3 = I^PP^In = N. 

Therefore, the transfer function from w to y is NC\(sI — A^^Bi. This completes the proof. □ 
Corollary 3 For any outer coupling matrix T, if the sum of every its row is 0, or the sum of every 

its column is 0, then the H2 norm of network ([3|) is 72-/V when it is stable, where 72 is the H2 norm 

of the single-node system (fJJ) . □ 
Remark 3 By Theorem 2 or Corollary 3, the H2 norm of system ^ with Yd given in ([SD or with 

any diffusive matrix is equal to 20 when N = 20 in Examples 1 and 2 discussed above, verifying the 

linear growth speed of the network energy. 

Next, consider networks with an antisymmetrical coupling matrix, 

= (jij)NxN, if* > j, 7ij = -1; if* < h lij = 1; and7ij = 0, i, j = 1, ■ • • , N. (10) 

Obviously, all the eigenvalues of Ta are located on the imaginary axis. Table 2 shows the H2 norms 
of network ([3]) with the outer coupling matrix V = Ya given by (jlOp . 



Table 2 H2 norms of network ([3]) with T — Ta for system data given in Section 3. 



N 


1 


2 


3 


4 


5 




10 




20 




Ex 1 


1 


1.8397 


2.5542 


3.1978 


3.8013 




6.6168 




12.1951 




Ex 2 


1 


1.7949 


2.4851 


3.1381 


3.7805 




7.0102 




13.5857 
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Similarly to the diffusive coupling case, one can see from Table 2 that the H2 norm of network Q 
with an antisymmetrical coupling also increases slowly. 

5 Networks with block-diagonal-input and block-diagonal-output ma- 
trices 

In this section, consider network ([3]) with block-diagonal-input and block-diagonal-output matrices, 
i.e., 

B = diag(.Bi, • • • , Bi), C = diag(Ci, • • • ,Ci), (11) 
but A remains unchanged. Table 3 shows the H2 norm changes of network §3§ with (|11|) . 



Table 3 H2 norms of network §5§ with (fTTj) and T given in (HJ) for system data given in Section 3. 



N 


1 


2 


3 


4 


5 




10 




20 




Ex 1 


1 


1.8974 


2.7203 


6.9089 


11.4222 




1.0026 x 10 3 




3.7236 x 10 6 




Ex 2 


1 


2.5531 


4.3147 


23.5579 


43.9595 




3.4090 x 10 4 




5.5723 x 10 9 





Comparing Table 3 with Table 1, one can see that the norms in Table 3 also increase dramatically. 
For small N, the norms in Table 3 are smaller than the ones in Table 1. However, they become larger 
than the ones in Table 1 when N is large. 



Table 4 H2 norms of network with (fTTj) and given in ([5]) for the same system data. 



N 


1 


2 


3 


4 


5 




10 




20 




Ex 1 


1 


1.1952 


1.4918 


1.5780 


1.7967 




2.2347 




2.9020 




Ex 2 


1 


1.1602 


1.4915 


1.5494 


1.8108 




2.2498 




2.9944 





Table 5 H2 norms of network Q with ((11)) and Ta given as in (| 10[) for the same system data. 



N 


1 


2 


3 


4 


5 




10 




20 




Ex 1 


1 


1.4142 


1.7321 


2 


2.2361 




3.1623 




4.4721 




Ex 2 


1 


1.4591 


1.8108 


2.1012 


2.3533 




3.3311 




4.7108 





Comparing Tables 4 and 5 with Theorem 2 and Table 2, one can see that the H2 norms of network 
([3]) with (jlip also increase slowly for the couplings and Ta, even much slower than the norms of 
the networks with column-input and row-output matrices. 

Corollary 4 For first-order node systems, the H2 norm of its corresponding network with anti- 
symmetrical coupling matrix Ta is 72V^V : . 
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Proof By Lemma 1, in this case, Q can be taken as a diagonal matrix in the form of al, and 
B T QB = j 2 N, where 72 is the H 2 norm of the single-node system. This completes the proof. □ 
Remark 4 One can clearly verify the result of Corollary 4 by examining Table 5. 



6 Influence of energy on nonlinear Lur'e networks 

The absolute stability of Lur'e systems and synchronization of Lur'e networks have been extensively 
studied (see 0[15] and references therein). In this section, consider the influence of the energy of the 
locally linerized system of a nonlinear Lur'e network on its network output. Given a Lur'e system as 
follows: 

± 1 = (A 1 -2A 12 )x 1 + B 01 f 1 (y 1 ), n(0) = B O i, 

(12) 

yi = Cqixi, 

where x\ is the state, y\ is the measured output, xi(0) is the initial condition, 

Ai = I ° ' I , B 01 = I ° I , C01 = ( 2 2 ), Au = B 01 C i, 



1 > 














1 C01 


-4 -2.5 1 




V 1 J 





referring to Example 2 in Section 3 for these matrix data, and the nonlinear function fi(yi) = \yi + 
1| — 1 2/1 — l|i which is just a piece wise linear function as used in the canonical Chua's circuit |15j . 
Obviously, fi satisfies the following sector condition: 

<M/i)< 2) / 1 ( ) = o. (13) 

yi 

Consider a network with regulated output formed by the node equation (|12p . as follows: 

f x = (I N <S)(Ai - 2Aia) + T L (g) A 12 )x + Bf(y), x(0) = E N <g> xi(0), 
y = Cix, (14) 
z = C 2 x, 

where x = (xf,--- ,arjy) T , X{ is the state of the i-th Lur'e system, y = (yi, • • • ,yv) T , /(y) = 
{fiiyi),'" Ji(yN)) T , B = I N ®B i, Ci = I N <g)C 01 , C 2 = £^®C i, 2 is the regulated output, 
-Eat is given as in ([3]), and is the outer coupling matrix similarly to T in ([3]). 

Linearizing network (|14p near the zero equilibrium, one gets a linear network as follows: 

x = (I N ®Ai + T L ®A 12 )x, x(Q) = E n ®x 1 (Q), 

(15) 

z = C 2 x. 

By the LQR control method [10], network (|15p can be viewed in another form as follows: 
x=(I N ®Ax + T L ®A 12 )x + E N ®x x {p)5{t), x(0-) = 0, 

(16) 

z = C 2 x, 



where 6(t) is the impulse function. It is well-known that the output energy of network (|15p is equal 
to the H2 norm of network (|16p. i.e., 

\\C 2 (sI - I N (g)A 1 - T L (g) A 12 )- l E N (g) xi(0)|| 2 . 

According to the discussions in Sections 2 and 3, with the matrix data in (|12p . the output energy of 
(|15p increases exponentially fast with = T as in @ and increases linearly with = as in 
([8]). Fig. 1 and Fig. 2 show the influence of energy changes of network (|15p or (|16p on the regulated 
output of Lur'e network (|14p . 




(a) The output of Lur'e (b) The output of network (fill) (c) The output of network (fl"4"|) 

system (fl"2l). with T L = T as in fll). with T L = r D as in 

Fig. 1 The output of a single-node system (I12p and outputs of Lur'e networks for iV = 16. 




(a) The output error of network (|14p (b) The output error of network (|14D 

with = T as in . with = To as in ([8]) . 

Fig. 2 The output errors of Lur'e networks with two initial conditions: 
x(0) = E N (g) cci(O) and x(0) = 0.95^jv (g) xi(0), for iV = 16. 
From Fig. 1 and Fig. 2, one can see that the shape of the solution of the diffusive network is 
similar to the solution of a single-node Lur'e system, but the output error of the diffusive network is 
smaller, for the two initial conditions shown in Fig. 2 (b). However, for the newly constructed coupling 
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relationship = T as in the network output is very different with the output of a single- node 
Lur'e system and the network output error with two initial conditions is larger. 

Remark 5 In fact, by the above idea, an input to output performance index can be established 
for network synchronization problems, especially for local synchronization problems. It is imaginable 
that the synchronization performance is bad for a nonlinear network when the energy of its locally 
linearized network is large. Therefore, large energy changes should be avoided for synchronization 
problems. 

7 Conclusion 

In this paper, the ^-norm energy accumulation problem has been addressed for linearly coupled 
dynamical networks. Three types of networks, i.e., a newly constructed network, a typical diffusively 
coupled network and an antisymmetrically coupled network, have been studied on the changes of their 
H2 norms with respect to the network size N. For the newly constructed network, the H2 norm 
increases exponentially fast, even much faster than 2 N when N is large. However, the H2 norms of the 
diffusively coupled and antisymmetrically coupled networks increase only linearly or even slower. This 
shows the complexity in energy changes of dynamical networks, despite the linear dynamical nature 
of the networks, which may be unexpected or even surprising. It has also been shown that large 
energy changes can have bad influence on network output, therefore should be avoided in network 
synchronization problems. Finally, it should be noted that the corresponding norm problem [5] 
can be similarly studied, which will be reported elsewhere. 
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